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1. Introduction 

Since tlic realization in experiments, Bose Einstein eondensate (BEC) of ultra-cold atoms 
has attracted much interest of theoretical physicists. Under optical dipole traps, the hypcrfine 
spin of atoms remains to be active and leads to the spinor BEC.^"^ So far, spinor BECs 
have been found to show a variety of phases. '^"^'^ For the hyperfine spin F = 1 state, the 
ground-state phase can be either polar (antiferroniagnetic) with ^"^Na^'^ or ferromagnetic 
with ^^Rb.^"^ F = 2 condensates have been realized with both ^^Na^ and ^''Rb^^^ atomic 
species. The F = 2 BECs arc classified into three distinct phases referred to as ferromagnetic, 
polar, and "cyclic". ""^^'^^ However, the ground-state phase of the F = 2 state at zero magnetic 
field is under discussion-*^^ due to its very short lifetime (a few milliseconds) for which the 
equilibrium state cannot be reached. Compared to the F = \ cases, the spin dynamics of the 
F = 2 BEC is less well-understood, especially in the cyclic phase, giving rise to experimental 
and theoretical challenges. 

Recently, solitons of spinor BECs in one-dimension have been studied analytically^^"^'' and 
numerically. ^^"^'^ In experiments, matter-wave dark and bright solitons are produced only for 
single-component BEC.^'^"^'' For a generic hyperfine spin F, the dynamics is described by 
the (2F -|- l)-component Gross-Pitaevskii (CP) equation. The one-component CP equation is 
called the nonlinear Schrodinger equation (NLSE). If all the spin-dependent interactions vanish 
and only the intensity interaction exists, the multi-component CP equation is equivalent to 
the vector NLSE, which is also called the Manakov equation. The soliton solutions of these 
systems are well-known. ^^"^-"^ For F = 1, at special sets of coupling constants one-soliton 
solutions and two-soliton collisions were explicitly shown for the bright soliton under the 
vanishing boundary conditions, ^^'^'^ and for the dark soliton^^ and the bright soliton^*^ under 
the nonvanishing boundary conditions, by finding the map from the GP equation to the 2x2 
matrix NLSE, which can be solved by the inverse scattering method. ^^'^^ However, for higher 
spins, such a map to a known integrable equation has not been found and we need to look for 
alternative methods. 

In this paper, we aim at seeking higher-spinor BEC solitons in one-dimension. For this 
purpose, we employ two different methods for demonstrating F = 2 BEC bright one-solitons. 
One facile method, the single-mode analysis, effectively utilizes the reduction of the multi- 
component GP equation to the one-component one. Another new method is a generalization 
of Hirota's direct method^^ to multi-components. Hirota's direct method has been successfully 
applied to get solitons in one-component equation systems, and we prove its applicability and 
strength even for the multi-component systems. Indeed, this method gives solitons beyond 
the single-mode analysis. Both methods are regarded to be effective even for nonintegrable 
equations, and suitable for the first step of investigation. The results include not only ordinary 
single-peak solitons but also twin-peak solitons which cannot be expressed as the superposition 
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of two singlc-pcak solitons. 

The paper is organized as follows. In Sec. 2 the GP equation for the F = 2 spinor BEG 
is introduced. In Sec. 3 we study the single-mode analysis. In Sec. 4 Hirota's direct method 
is generalized for the F = 2 GP equation. In Sec. 5 we present one-soliton solutions obtained 
by Hirota's method and discuss their properties. The last section is devoted to discussion and 
conclusion. 

2. F — 2 Spinor BEC in One-Dimension 

In the mean-field theory, the F = 2 spinor BEC is characterized by the local order param- 
eter (or, the macroscopic wavefunction) with five components, * = ($2,^1,^0)^-1)^-2), 
reflecting the five spin degrees of freedom. For the magnetic quantum number j = —2, - ■ ■ ,2 
with respect to the quantization axis chosen in the 2;-direction, $j = ^j{x,t) = {^j{x,t)). In 
words, $j are given by the ground state expectation value of the boson operators 4/j{x,t), 
which satisfy the equal-time commutation relation [^aix,t),'^j^{x' ,t)] = SapS{x — x') for 
a,/3 = -2,--- ,2. 

We consider the dynamics of the F = 2 spinor BEC in one-dimension. The evolution 
equation for the local order parameters is described by the multi-component Gross-Pitaevskii 
(GP) equation, 

dt (5** ^ ' 

Here the energy functional is defined by^"^^ 

The coupling constants Cj are real and can be expressed in terms of a transverse confinement 
radius and a linear combination of the s-wave scattering lengths of atoms. The interaction 
energy is derived from the short-range interactions of atoms in the scattering channel with 
total spin 0, 2, 4, and is given in terms of the number density 



n 



the spin densities / = (/^, f^, /^), where for i = x,y, z, 



r= E ^«fa/3^/3' (4) 
a,/3=-2,---,2 

and the singlet-pair amplitudc^^'^^ 

e = 2$2^-2 - 2$i$_i + (5) 

The meaning of 9 is clear if we write with the Clebsch-Gordan coefficient as 6 = 
V^Sj j'(00|2j; 2_7')$j$jv, i.e. it measures the formation of spin-singlet "pairs" of bosons. 
The prefactor \/5 is introduced just for convenience. The C4-term in (2) includes the scatter- 
ing process 2 -|- (—2) -(-^^ -|- 0, which changes the z-component of the spin states of bosons 
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by two and is absent for F = 1. We also write the spin densities as 
matrices f* in = 2 are explicitly represented as 
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Without the magnetic field, the energy (2) is invariant under an SU(2) rotation and the system 
has an SU(2) symmetry. In particular, an obvious symmetry is the one under exp[i7rf^] : $j i— > 

We set h = 1, 2m = 1 to simplify the expressions. With an operator L = \dt + d"?;., the 
explicit form of GP equation (1) is 

L^±2 = con$±2 + C2 (±2/^$±2 + /^^±i) + C4G$;2> (7a) 



£.$±1 = con$±i + C2 I /±$±2 ± r^±i + ^/^$o I - C4e$;i, 



V6 



L<^o = con^o + C2 ( + ^/"^-i ) + C40^o- 



(7b) 



(7c) 



The right-hand sides include cubic terms with respect to $j . If the spin-dependent interactions 
are absent, i.e. 02 = 04 = 0, the GP equation is reduced to the Manakov equation with five 
components and solutions for the initial problem as well as multi-solitons are known in the 
formalism of the inverse scattering method.^°'^^ However, this is a trivial reduction for spinor 
condensates. In the presence of the spin-dependent interactions where rich phenomena are 
expected, the GP equation becomes highly correlated and hard to be solved explicitly. 

We concentrate on soliton solutions for the F = 2 GP equation. A soliton propagates 
keeping its own wave properties. Through its free translational motion, physical quantities 
given by the integral of densities characterize the soliton, such as the particle number N = 
J ndx, the spin F = J fdx and the volume of the singlet-pair S = J \&\dx. 

In the subsequent sections, we attempt to derive one-soliton solutions with non-trivial spin 
degrees of freedom and clarify their physical properties. We apply two methods and obtain 
several one-soliton solutions for the GP equation (7). 
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3. Single-Mode Analysis 

The single- mode analysis assumes the following amplitude for the order parameters: 

^{x,t) = Acl,{x,t), (8) 

where A = {A2, Ai,Ao, A^i,A-2)- That is, the order parameters have the same spatial profile 
but can have different magnitude. Normalization is such that Ylj l^jP = 1- We require the 
GP equation to lead to the one-component nonlinear Schrodinger equation for 0, 

idtcp + dl4> - c\(p\'^4> = (9) 

with C being a real constant. This imposes the consistency conditions on the nonlinear terms 
of the GP equation. By the freedom of SU(2) rotation, we can fix the spin in the z-direction 
to have = f~ = 0. Then the conditions read 

En = --- = Ej„ (10) 

for Ji, • • • , jfe G {-2, • • • ,2} with Aj^,--- , Aj,^ / 0, where 

^2 = 202/^ + 046^12/^2, (11a) 

El = - C4eA*_JAu (lib) 

Eo = c^eAl/Ao, (11c) 

E-i = -C2r - c^SAl/A^i, (lid) 

E_2 = -2c2r + C4@A*/A_2- (lie) 

Here, P = Y.a,p and 6 = 2A2A_2 - 2Ai^_i + Al for simplicity. 

We summarize the result of the examination of the consistency conditions. It is sufficient 
to specify amplitudes by their representatives by virtue of the SU(2) symmetry. 

• The ferromagnetic states, |/| > and |6| > 0. 

A= (p2,0,0,0,p_2), 4C2=C4. (12) 

A= (0,pi,0,|^_i,0), C2 = C4. (13) 

^ = (^'2,Pl,P0,P-l,P-2), C2 = 0,6 = 0. (14) 

• The polar states, f = and |0| > 0. 

A= {p2,0,po,0,P-2), \p2\ = \P-2\, C4 = 0. (15) 

A={p,q,r,-q*,p*), p,qeC,reR. (16) 

• The cyclic state, / = and |0| =0. 

A={p,0,0,V2q,0), p,qeC, \p\ = \q\. (17) 
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Here wc can take pi G C. These three spin states for F = 2 spinor BEC are specified in Ref. 
11. The cyclic state, which is absent in the F = 1 case and is available for spin F > 2 boson 
systems, exhibits unusual features such as phase-locking phenomena and kink excitations^^ 
owing to a unique nature that the condensate energy depends on the relative value among the 
phase factors of ^j. It is remarkable that any onc-solitons in the three states are obtained. 
The coefficient of the nonlinear term in (9) turns out to be 

C = Co, for e = 0, (18) 

C = cq + C4, otherwise. (19) 

When the effective coupling is the attractive one, C < 0, we have the bright one-soliton 

4>{x,t) = ^e'^^sechxr. (20) 



The position function Xr and the phase function Xi of the soliton are given by 

Xr = "^krkit — kiX + 6r, (21) 

Xi = -{kl - ki)t + krx + Si, (22) 

respectively, where x = Xr + ^Xi = kx — k'^t + 6 with k = kr + iki. One can also see the plane- 
wave solution (t){x,t) = exp[i(ii'a; — fit)], where K and $7 are real with $7 = K'^ + C. Under 
C < 0, this plane-wave is unstable against the modulation and is decomposed into bright 
solitons during time-evolution. In the case of the repulsive coupling C > 0, a dark-soliton is 
formed under the nonvanishing boundary conditions \(f)\ — ^ const, as a; — ^ ±oo. 

We remark on the reduction to F = 1. If C4 = and * = (0, $o/\/3, 0) such that 
+ $Q = 0, the GP equation for is equivalent to the F = 1 GP equation. The 
polar single-mode soliton with (16) and p = is reduced to the F = 1 polar one-soliton in 
Refs. 16-18. 

4. Hirota's Direct Method 

In this section, we introduce Hirota's direct method. Hirota's direct method is powerful 
for getting solitons in both integrable and nonintegrable one-component partial differential 
equation (PDE) systems. We generalize this method for the F = 2 GP equation. 

By putting $j = Gj/H for j = —2, • • • ,2, the GP equation is transformed into the form, 

(.A + i.aO..«-(||J + |§|)«^^0, ,23a) 
dIH-H + cq ^ |G„|^ = 0. (23b) 

a=-2,-,2 

Here the Hirota derivative is defined as 



dti dt2 1 \ dxi dx\ 



V2 



(24) 

*1— *2— * 
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In the Manakov and the one-component systems, the C2- and C4-terms in eqs. (23a) are absent 
and the transformation is called the bilinear transformation. 

For attractive spin-independent interaction, we can set cq = —2 by scaling the order 
parameters. Equations (23) give solitons with the vanishing boundary conditions \^\ ^ as 
X zbco, that is, bright solitons. For the nonvanishing boundary conditions, extra terms are 
needed. 

Within Hirota's method, a one-soliton solution may be obtained by a finite-order pertur- 
bation, 

H = l + e'^h2 + e'^h4, (25a) 

Gj = sgij+e^g3j. (25b) 

Substituting (25) into (23), we solve them order by order. At the first order, we have 

gij = n^e'^*^^-*^'*) (26) 

with k = kr + iki e C and Uj {j = —2, • • • ,2) being free parameters for one-solitons. A 
success of terminating the perturbation expansion at a finite order leads to soliton solutions. 
We cannot expect soliton solutions for generic values of parameters. Our strategy is that by 
utilizing the freedom of 7 parameters, i.e. 5 for one-soliton Uj and 2 for interaction couplings 
C2 and C4, we look for both one-soliton solutions and the valid parameters in order to have 
those solutions. 

It is known that for one-component PDEs, Hirota's direct method is applicable for upto 
two-solitons even in nonintegrable systems.^^ In fact, the expansion (25) is usually taken for 
a two-soliton. In our case, as we will show later, we have a twin-peak one-soliton for the 
F = 2 GP equation. In a sense, a twin-peak one-soliton may be considered as a degenerate 
two-soliton. This kind of observation is also seen in Refs. 19, 33. For a general two-soliton, we 
take (26) as a combination of two plane waves. In principle, those calculations which include 
more higher orders and multi-solitons are possible in integrable systems, but we do not reach 
the argument of integrability and postpone them to future works. 

5. Results for One-Solitons 

We present one-soliton solutions through Hirota's direct method. A one-soliton is deter- 
mined by the following parameters; fc^: (half of) the phase velocity of the envelope soliton, 
kf. the amplitude of the soliton, 11 = (H2, Hi, Hq, H_i, n_2): distribution among spin com- 
ponents. The position function and the phase function of the soliton are the same as the 
one-component one's, (21) and (22), respectively. We normalize 11 by = 1 since its 

factor can be absorbed in the shift of x- 
5.1 Single-mode (Single-peak) soliton 

Single-mode one-solitons are reproduced with specific values of coupling constants. 
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a) A single-mode soliton in the ferromagnetic state, with |/| > and = 0. 

n = {p2,pi,po,p-i,p-2), Pi e c, 2p2P-2 - 2pip^i+pl = 0. 
C2 = 0, C4 = arbitrary. 

b) A single-mode soliton in the polar state, with / = and |0| > 0. 

^= iP2,0,P0,0,P-2), Pi&C, \p2\ = \p-2\- 
n={p,q,r,-q*,p*), p,qeC,reR. 
C2 = arbitrary, C4 = 0. 

c) A single-mode soliton in the cyclic state, with / = and = 0. 

n= (p,0,0,V2g,0), p,qeC, \p\ = \q\. 
C2 = arbitrary, C4 = arbitrary. 
For the cases a)~c), the one-soliton solution is 

* = fcie^-^'sechxrll. 

It should be remarked that the cyclic soliton in c) exists for all interaction couplings. 

d) A single-mode soliton in the polar state, with / = 0, |0| > and a doubled particle number 
compared to the cases a)^c) for the same kr and ki. 

n = (p,g,r,-g*,p*), p,qeC,re'R, 

n = (p2,o,po,o,p-2), Pi G c, \p2\ = \p-2\, pIIp2P-2 e M. (3i) 

C2 = arbitrary, C4 = 1. 

e) A single-mode soliton in the ferromagnetic state, with |/| > 0, |0| > and a doubled 
particle number compared to the cases a)^c) for the 



(29) 



(30) 



n = (p2,o,o,o,i?_2), Pi^ 

C2 = 1/4, C4 = 1. 



(32) 



n = (p,o,o,o,o), pec. 

C2 = 1/4, C4 = arbitrary. 



(33) 



n = (o,pi,o,p_i,o), kgc. 

C2 = 1, C4 = 1. 



(34) 



n = (o,p, 0,0,0), peC. 

C2 = 1, C4 = arbitrary. 



(35) 
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For the cases d) and e), the one-sohton solution is 

* = V2kie'^'sechxr'n. (36) 

Compared to (30) in the previous cases, the amplitude is larger by a factor \/2. Prom the 
single-mode analysis in Sec. 3, it is shown that the effective coupling C for the one-component 
equation has the value cq + 04 = —1, not cq = — 2. Accordingly, the amplitude gets the factor 
\/2, and hence the particle number becomes doubled. 

The single-mode solitons in this section are included in the result in Sec. 3. 

One can see that for single-mode one-solitons, zero local spin is allowed for arbitrary C2, 
and zero singlet-pair amplitude is allowed for arbitrary C4, since the corresponding interaction 
energy in (2) is ineffective, respectively. Moreover, the solitons in the cases a)~c) turn to have 
vanishing interaction terms with C2 and C4 in the energy (2), and they are regarded as the 
solitons in the Manakov system. 
5.2 Twin-peak soliton 

We further investigate one-solitons which cannot be expressed within a single-mode form. 
They have a wave-form with twin peaks. The distance of the twin peaks is freely adjusted by 
changing the parameters of the one-soliton. One may be tempted to say that the twin-peak 
soliton is the superposition of two identical single-peak solitons with shift of their positions, 
but it is not true because physical densities of the twin-peak soliton are not always just the 
sum of those of two single-peak solitons. Such a twin-peak one-soliton was already discovered 
in the F = 1 spinor BEC,^^'^^ but is allowed only for the polar state, i.e. with zero total spin. 
In our result for F = 2, we find that twin-peak one-soliton occurs both in the polar state and 
in the ferromagnetic state. 

/) A twin-peak soliton in the polar state, with F = (but locally / 7^ 0) and |B| > 0. 

n = (p2,pi,po,p-i,p-2), Vi e c. ^^^^ 

C2 = 0, C4 = 1. 

The wavefunctions have the form 

V2fc, (-l)VpV^'-+Pje- >^'' 
cosh 2xr + cosh u) 

for j = —2, • • • ,2, where T = 2p2P-2 — "^PiP-i + Pq, = cr = ^/|^| S'ld coshw = 

N/ITI. Physical densities are calculated as follows; 



= ^ 17^ , IT., e'^^ > (38) 



n = kf 



sech^ 



+sech2(xr + |)] , (39) 



^ * |T| (cosh 2xr + cosh 6c;^2' 



4A-2 

lei = , (41) 

cosh 2xr + cosh uj ' 
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Fig. 1. Density plots for /) (twin-peak polar soliton) with 11 = (1, 3, 3, 1, 3), fc^ — h ~ 1. (a) Densities 
for each component (|<i>2p: red, purple, |<I>op: blue, yellow, |<i>-2p: green), (b) The 

number density, (c) The spin densities (/^: green, blue, pink), (d) The absolute value of 
the singlet-pair amplitude. 



where f = X^q /3 n* f^^II/j. The total amounts of the quantities are obtained by integrating 
these densities as 

N = 4h, (42) 
F = (0,0,0), (43) 
S = AkiU) cosechu;. (44) 

Figure 1 shows an example of the twin-peak polar soliton. We observe that the spins 
contained in the two peaks have the same amount with the opposite sign, and form a polar- 
ization. Therefore, in total the soliton has zero spin. The reduction to a single-peak soliton is 
achieved by sending 'J ^ 0. In this limit, the two peaks get infinitely far apart and eventually, 
the remained single-peak soliton contains nonzero total spin and no singlet-pair amplitude, 
and coincides with the ferromagnetic one (14) in a). We also observe that the singlet-pair 
amplitude is localized around the center of the twin peaks. This indicates that the twin-peak 
soliton is not just the superposition of two identical single-peak solitons. The same argument 
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also holds for the next case g). 

For the special case with zero local spin, where we make f = for i = x,y, z, the following 
set of parameters is allowed: 

n = (p2,0,po,0,p-2), Pi&C, |p2| = |p-2|- 
C2 = arbitrary, C4 = 1. 



(45) 



The wavefunctions are 



P±2 



(2p2P-2 + Po)^ + cosh LO 

$±1 = 0, 



$n = 



(2p2P-2 + Po)^ cosh 2xr + cosh a; 



where cosh a; 



2|p2P-2| + b_0p 
|2p2P^2+Pol 



^ = arg 1 + 



2p2P-2 



and ^' = arg ( 1 + 



2p2P-2 



(46a) 
(46b) 
(46c) 
In particular, 



$, = ^'{= 0) (coshu; = 1) gives the single-mode soliton of d). 
g) A twin-peak soliton in the ferromagnetic state, with |/| > and \Q\ > 0. 

n = (p+,o,o,o,p_), pi^c. 

C2 = 1/4, C4 = 0. 

n = (o,p+,o,p_,o), piGC. 

C2 = 1, C4 = 0. 

We use s = 2 for the case (47) and s = 1 for the case (48). Then, the wavefunctions are 

- |;>_|2| cosh2xr + coshw 

and the others are constantly zero, where a = sign (|p+P — and cosh a; = ijp^pl'^jp'pi 

The densities are given by 



(47) 



(48) 



(49) 



n = k7 



sech' 



-I- Sech"' Xr + TT 



r = p = 0, 

4fc?SC7 

/ 



lei 



cosh 2xr + cosh uj 
Akf I sinh u; sinh 2xr 



(cosh 2xr + coshti;)^ 
The total amounts of the quantities are 

N = Aki, 

F = {0,0, Akisaoj cosechw), 

P- 



S = 4L- 



P+ 



(50) 
(51) 
(52) 

(53) 

(54) 
(55) 

(56) 
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(c) (d) 

Fig. 2. Density plots for g) (twin-peak ferromagnetic soliton) with 11 = (0, V5/2, 0, 1, 0), kr ^ ki = 1. 
(a) Densities for each component purple, yellow), (b) The number density, (c) The 

spin densities (/^: green, f^: blue, pink), (d) The absolute value of the singlet-pair amplitude. 



Figure 2 illustrates an example of the twin-peak ferromagnetic soliton. One can see that 
the spin densities are localized around the center of the twin peaks. As — > 1, the two 
peaks of the soliton get infinitely far apart and a single-peak soliton in the polar state with 
(15) in b) is left. It is interesting that the reduction of the twin-peak soliton to the single-peak 
one changes the state from ferromagnetic to polar in g) and vice versa in /). 

Note that (49) is the solution of the two-component coupled NLSE,^^ 

= -(a|$,|2 + /3|$_,|2)$„ (57a) 

= -(/3|$,|2 + a|$_,|2)$_„ (57b) 

with a = 1 and /3 = 3. It was established that two-component coupled NLSE is integrable 
only for a = /3, corresponding to the original Manakov equation.'^'' 

6. Discussion and Conclusion 

We have studied one-soliton solutions for the Gross-Pitaevskii (GP) equation of the F = 2 
spinor Bose-Einstein condensates (BEG) by means of two methods, the single-mode analysis 
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and the multi-component generalization of Hirota's direct method. The latter method has 
been successfully applied to show twin-peak solitons both in the ferromagnetic and the polar 
states, which cannot be accessed by the single-mode analysis. 

Hirota's method is not restricted to the present analysis. One can also find soliton solutions 
for higher-spinor BECs or other types of multi-component systems by generalizing this method 
as presented in this work. For instance, applying to the F = 1 spinor BEC, we reproduce the 
bright solitons in Refs. 16-18, at the order e"^ for the single-mode solitons and at the order 
for the twin-peak polar soliton. Mathematically, incorporating such as C2- and C4-terms of 
eqs. (23a) in Hirota's method suggests a new direction for further extension of the framework. 

One of our next interests is the integrability of the GP equation. In the integrable systems 
with multi-components, multi-solitons should be constructed from any combinations of one- 
solitons. Their collisions are factorized into successive two-soliton collisions, but in contrast to 
one-component systems, it is not always the case that each soliton keeps its shape after colli- 
sions, deforming its parameters for internal degrees of freedom. The results in this paper pick 
up specific interactions, and we hope that they especially include integrable points. Whether 
those systems are integrable or equivalent to already known systems is an interesting future 
problem. The Painleve analysis may give a clue for the problem. 

From the physical point of view, the discovery of all one-solitons in the ferromagnetic, polar 
and cyclic states is of much importance. As higher-spinor BECs should exhibit richer physics, 
we expect that higher-spinor solitons will make wider possibilities in various applications. 
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